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ABSTRACT 

During the past decade or so, measurements of Galactic HI absorption using VLBI 
against extra-galactic sources, as well as multi-epoch observatios in pulsar directions, 
have detected small-scale transverse variations corresponding to tens of AU at the 
distance of the absorbing matter. Hitherto these measurements have been interpreted 
as small-scale structure in the HI distribution with densities uhi ~ 10 4 — 10 5 cm~ 3 , 
orders of magnitude greater than those of the parsec-scale structure. Naturally it is dif- 
ficult to imagine how such structures could exist in equilibrium with other components 
of the ISM. 

In this paper we show that structure on all scales contributes to the differences 
on neighbouring lines of sight, and that the observed differences can be accounted for 
by a natural extension of the distribution of irregularities in the distribution of HI 
opacities at larger scales, using a single power law. This, in our opinion, should put an 
end to the decades long puzzle of the so-called small-scale structure in HI and other 
species in the Galaxy. 

Key words: interstellar medium: clouds — interstellar medium: structure — radio 
lines: atomic — interstellar medium: molecules — pulsars: general — interferometry: 
interstellar radio sources: 21 cm radiation 



1 INTRODUCTION 

The warm component of the Galactic neutral atomic hydro- 
gen (HI), studied extensively through its 21-cm line emission 
using single-dish measurements, shows largely uniform dis- 
tribution and has revealed the large-scale structure of our 
galaxy. The cold atomic component has also been probed us- 
ing the 21-cm absorption observable in the spectra of bright 
continuum sources in the background. The earliest interfer- 
ometric study towards the bright supernova remnant Cas-A 
(Clark, Radhakrishnan & Wilson 1962; Clark 1965) revealed 
the presence of structure finer than known earlier in cold HI. 
Later, many aperture-synthesis observations of the Perseus- 
arm features detected structures down to the resolution limit 
(an arc-minute) in these observations (Greisen 1973; Breg- 
man et al. 1983; Schwarz et al. 1986). These and various 
other indications had suggested that the HI gas in our galaxy 
is organized on, and maintains, a hierarchy of scales from 1 
kpc to at least 1 pc. Contribution from scales much smaller 
than the parsec-scale was expected to be a tiny fraction of 
the total (see, for example, Dickey & Lockman 1990). There 
were no serious difficulties in understanding this picture. 

However, subsequent HI absorption studies using 
VLBI observations of extra-galactic sources (as background 
sources) triggered what has remained as a puzzle for a few 



decades. Dieter, Welch & Romney (1976) were the first to 
note variations in the HI opacity on a scale smaller than 
0.16" (in the direction of 3C147). This they interpreted as 
implying a cold-HI cloud size smaller than 70 AU and a 
volume density in the cloud of ~ 10 s cm -3 . More than a 
decade later, Diamond et al. (1989) reported more VLBI 
observations on 3C147 as well as two other extra-galactic 
sources, supporting the conclusions of Dieter, Welch & Rom- 
ney (1976). They interpreted their findings as suggesting 
linear diameters of the absorbing clouds to be as small as 
~25 AU and correspondingly high densities. More elaborate, 
recent VLBI observations (Davis, Diamond & Goss 1996; 
Faison et al. 1998; Faison 1999) confirm some of the ear- 
lier reports of opacity variations on small transverse scales, 
while in some other cases find no detectable variation across 
different components of the background sources. 

A new technique employing multi-epoch Hl-absorption 
measurements in pulsar directions was suggested by Frail et 
al. (1991) as well as by Deshpande et al. (1992) indepen- 
dently. Based on an extensive multi-epoch study of seven 
pulsars (sampling various directions in the galaxy), Frail et 
al. (1994) reported optical-depth changes of <0.1 on the 
small spatial-scales of 5 AU to 100 AU and concluded that a 
significant fraction (10%- 15%) of the cold HI gas is in such 
small-scale structure. 
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Naturally, there are serious difficulties about such a 
structure being in pressure equilibrium with the other com- 
ponents of the medium given the estimated volume den- 
sity being so high, and hence about what processes would 
generate and help maintain such apparently commonly en- 
countered structure. In a recent paper, Heiles (1997) has 
summarized the main results from these observations of the 
so-called small-scale structure in the interstellar HI (as well 
as from optical observations of interstellar absorption lines 
of Nal & Call). The paper also points out that existence of 
such a tiny-scale atomic structure would imply, under con- 
ventional interpretation, over-abundance of H2 leading to 
very large extinction. To ease such difficulties, Heiles (1997) 
has proposed geometric solutions invoking structures con- 
sisting of cold, dense curved filaments or sheets (that line-up 
along sight-lines) to explain the observed variations in the 
Hi-opacity, but with moderate values of the implied volume 
densities. On the other hand, Dickey & Lockman (1990), 
based on many arguments, conclude that while small-scale 
structure does exist, it is only a tiny fraction of the total HI 
column density along any sight-lines. 

The two important inter-related questions raised by the 
apparently puzzling detections of the opacity variations over 
small spatial scales, and often asked, are the following, a) 
Does the atomic medium resemble the diffuse ionized com- 
ponent ? That is, does it have a power-law distribution of 
sizes ? If so, is the behavior similar over the whole range of 
scales probed ? b) Is the AU-sized structure only peripher- 
ally related to the parsec or larger scale structure in HI ? 
That is, does it represent a physically distinct population of 
structures ? 

In this paper we show that the observations have been 
misinterpreted, and that the observed small-scale structure 
is not at all unexpected. The observed opacity differences are 
consistent with a single power-law description of the distri- 
bution of HI opacities in the interstellar medium, and almost 
all scales contribute significantly to the observed differences. 
It was incorrect to assume that the observed structure must 
be due to high-density clouds whose longitudinal dimension 
is the same as the separation between the lines of sight at 
distances comaprable with the distances of the absorbing 
matter. 



2 WHAT DO WE ACTUALLY MEASURE IN 
THE VLBI AND MULTI-EPOCH PULSAR 
OBSERVATIONS ? 

2.1 Effects of the source structure and telescope 
filter function 

The situation in the Hi-opacity-variation measurements us- 
ing VLBI on extra-galactic sources and the multi-epoch 
pulsar observations^] can be analyzed by identifying three 

* For further discussion, we will treat, without loss of any details, 
the multi-epoch pulsar measurements as equivalent to those made 
against a background source consisting of as many number of 
incoherent component sources as the number of epochs, and each 
component location defined by the apparent pulsar-direction at 
the corresponding epoch. Due to the interstellar scattering, the 
size of the component sources may appear larger compared to that 
intrinsic to the pulsar radiation. Although the pulsar observations 



important ingredients that dictate what actually we would 
measure; namely, 1) structure of the background source, 2) 
structure in the absorbing gas, and 3) the spatial frequency 
filter function of the telescope. As for their frequency depen- 
dence (over the observing bandwidths that are usually small 
compared to the centre frequency) , the first quantity can be 
assumed to be constant, the third one will vary only slightly 
but predictably, while the second quantity can vary consid- 
erably from channel-to-channel in frequency (or velocity). 

In the continuum channels of the observed band (triv- 
ial case of zero-opacity), the source structure apparent to 
the observer is of course the "true" structure of the back- 
ground source. Whereas in the line channels, which are the 
ones of interest, the apparent source structure is modified 
by the opacity structure. The apparent structure is always 
a product of the two "true" structures, making it, in gen- 
eral, appear finer (and consequently extending its visibility 
range to higher spatial-frequencies) compared to the indi- 
vidual ones. Two instructive cases are; namely, 1) Uniform 
finite opacity and VLBI-scale structure in the background 
source and 2) small-scale structure in opacity and uniform 
brightness background source. And, let us view these two 
situations, say, using one VLBI baseline. In the first case, 
even though the opacity is uniform (i.e. has no structure), 
the absorption will be visible as long as the "true" source 
structure (which it mimics) has finite visibility at a given 
baseline. Thus, the absorption "visibility" on a given base- 
line does not necessarily suggest a structure in opacity on 
the corresponding angular (or the related spatial) scale. In 
the second case, the continuum visibility would be zero, but 
there may be finite visibility in the line channels. So, any ob- 
served "visibility" should be directly attributed to a struc- 
ture in opacity. In fact, it provides a reasonably pure mea- 
sure of the power at the corresponding spatial scale in the 
spectrum characterizing the distribution of opacity (partic- 
ularly at small optical depths^). Green (1993) has indeed 
made such interferometric measurements of the 21-cm emis- 
sion line to directly sample the associated power spectrum 
at discrete spatial frequencies and Lazarian (1995) has pre- 
sented a technique to study the underlying 3-D character- 
istics using such measurements. We will return to this case 
again later. 

In reality the situation is somewhere in between the 
two relatively simple extremes and hence needs even more 
care in interpreting the measurements. An essential step 
then involves proper imaging of the apparent structures in 
continuum- and in line-channels and using the comparison 
to estimate opacity in the usual way. The opacity can be 
estimated only over the extent of the background source im- 
aged. The opacity estimates in closely spaced directions can 

are made usually with a single-dish, for the above equivalence to 
be complete, we consider the angular resolution to match the 
scatter broadened size. 

t Although we have used absorption and opacity as analogous 
to each other, what the measurements respond to is the modified 
source structure. And the fractional difference in the apparent 
structures in two spectral channels gives the structure in frac- 
tional absorption from which the opacity is to be computed. In 
this context, the particular non-linear correspondence between 
the depth of absorption and the opacity (optical-depth) should 
be borne in mind, particularly at large opacities. 
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be compared as is done in the studies using pulsar-probes 
and in more recent VLBI studies (Davis, Diamond & Goss 
1996; Faison et al. 1998; Faison 1999) that have carefully 
imaged the opacity distributions. 

2.2 Expected differences in the opacity & 
contributing scales 

Now, we ask and try to answer two crucial questions, a) 
What is the magnitude of opacity differences that we would 
expect to observe between a given pair of sight-lines ? and 
b) what scale(s) from the opacity distribution should be con- 
sidered as contributing to the observed opacity differences 
between two given sight-lines ? 

Let T v (x,y) represent a two-dimensional distribution of 
opacity in the transverse coordinates (x,y) for a given ve- 
locity channel. For simplicity, let us choose the transverse 
spatial separation corresponding to the angular separation 
between a given pair of thin sight-lines at the Hi-screen dis- 
tance to be along the x-axis and denoted by x . The rms 
value of the opacity difference (At v (x )) expected to be ob- 
served is given simply by the square-root of the structure 
function of r v (x,y) at a spatial scale of x . Analytically, 
this can be expressed as 

(At v (x )) 2 = S Tv (x ) = < (t v (x, y) - t v (x - x a , y)) 2 > 
where <> denotes ensemble average of the quantity over all 
(x,y). 

To examine the contributing scales to this opacity dif- 
ference at separation x , we consider the power spectrum 
(Pr(fx, fy)) as a function of the spatial frequencies (f x ,fy) 
(corresponding to spatial scales 1/ fx, associated with 

the distribution r v (x,y), such that P T is the Fourier trans- 
form of the autocorrelation function of r v (x,y). If the aver- 
age power spectrum can be described a s a powe r- law, i.e. 
< Pr(f 3 ) >= P?Js a (where /, = y/jf+J 2 and such 
that a is positive) then for 2 < a < 4, the structure 
function can also be described as a power law such that 
S Tv (x ) = S° v x a ~ 2 (Lee & Jokippii 1975). Then from the 

a — 2 

above equation, it follows that Ar v (x ) — At°x 2 . 

The important thing to note here is that the quantity 
At v (x ), which relates statistically to the opacity differences 
observers have measured, has a much slower dependence on 
the transverse separation (x ) between the sight-lines than 
the amplitude at the spatial frequency (\/x ) would have in 
the spectrum. Now, to make some quantitative estimates, 
we need to know the details of the power spectrum. Fortu- 
nately, such details are now available from a recent power 
spectrum analysis of opacity in the direction of Cas-A (Desh- 
pande, Dwarakanath & Goss 2000; hereafter DDG) using 
the opacity images measured by Bieging, Goss, & Wilcots 
(1991; hereafter BGW). The authors (DDG) report that the 
power spectrum is of a power-law nature over scales rang- 
ing from ~0.02 pc to ~4 pc and the value of a to be close 
to 2.75, significantly different from the Kolmogorov value of 
11/3. When viewed over 0.5 km/s wide velocity channels 
(similar to that used in the small-scale structure studies 
being discussed here), the rms variation in opacity across 
these images is about unity, making the structure function 
at x ~ 4 pc (corresponding to the angular size of Cas-A 
and the location of the absorbing cold HI) equal to about 
2. From this calibration and the value of a as suggested by 
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Figure 1. An azimuthally averaged power spectral version of a 
simulated 2-d spectrum of opacity distribution with a = 2.75 as 
suggested by the data in the Cas-A direction (DDG 2000). The 
y-axis is calibrated to indicate the power contributed by each 
spatial frequency to the r distribution. 

the Cas-A data, it follows that optical depth differences of 
typically 0.2 (rms) should be in fact expected at transverse 
separations of ~1000 AU. Given the power-law index of the 
structure function, at 100 AU separation, the expected (rms) 
differences would drop by a factor of ~2. To assess further 
the expected opacity differences, we should examine ideally 
the (probability density) distribution of the expected magni- 
tudes of differences for each value of x , the structure func- 
tion itself representing the second moments of such distri- 
butions as a function of x . Unfortunately, the observations 
of BGW have 7 smoothing which limits the range of scales 
(at the smaller-scale end) that we would like to examine. 
Hence, simulations avoiding such smoothing were consid- 
ered (see Fig. 1). A more detailed discussion on this and 
related issues is given in DDG (2000). A complex hermitian 
symmetric spectrum was simulated with the contributions 
(the real and imaginary parts) at different spatial frequency 
f s represented by uncorrelated random numbers following 
zero-mean Gaussian statistics having variance matched to 
the fg~ a power law. Such a spectrum over the 512x512 ma- 
trix (in the 2-d case, and a 2 18 point array in the 1-d case) 
was Fourier transformed to obtain a simulated r-distribution 
and suitably scaled to have an rms of unity. The spatial ex- 
tent of the distribution is assumed to correspond to ~ 4 
pc, consistent with the data in Cas-A direction. Using the 
simulated distribution of opacity (similar to that observed 
by BGW in the Cas-A direction), we have estimated, as a 
function of the spatial separation, two quantities indicative 
of the related (one-sided) probability distribution of the ab- 
solute differences. These are the rms value (i.e. the square- 
root of the structure function) and the maximum value of 
the magnitudes of opacity differences. The result of this es- 
timation is shown in Fig. 2 and it is clear that occasionally 
the opacity difference can be nearly an order of magnitude 
higher than the rms values (related to the structure func- 
tion). This makes the detected differences hardly surprising 
and therefore, they should be treated as only consistent with 
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Figure 2. Two quantities (related to the magnitude of differences 
expected as a function of the transverse scale probed) estimated 
from 1-d simulations of optical-depth (t) distribution are shown. 
The trends with symbols o,x indicate respectively the rms (i.e. 
square-root of the structure function) and the typical peak differ- 
ence magnitudes. The 1-d distribution of opacity over 5 orders of 
magnitude in scales (1:2 18 ) was obtained from a simulated spec- 
trum in 1-d with a power-law index of 1.75 (corresponding to a 
2-d equivalent index of 2.75). Due to the computational require- 
ments being too high, the 2-d case was not attempted on the 
same scale-range. However, the 2-d simulation referred to in Fig. 
1 were used to confirm consistency of the results over a range 
of longer scales spanning two orders of magnitude. The Y-axis 
values here are calibrated such that the maximum scale may be 
equated to 4 pc (see text for details). The r distribution used 
was obtained by Fourier transforming a simulated spectrum in 
the spatial frequency domain; The repetitive nature implicit in 
the Fourier transforms can affect estimation on scales close to the 
transform length. Hence, the above estimates are made only for 
scales that are less than half of the transform length. 
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Figure 3. A sample section from our (1-d) simulated opacity 
distribution is shown. The axis scales are calibrated consistent 
with the observations in the Cas-A direction (corresponding to 
the scale range at longer-scales). 



a single power-law spectral description of the opacity distri- 
bution (e.g. as derived in the Cas-A direction). Fig. 3 shows 
a small section of the simulated distribution as an exam- 
ple of the expected variation on small transverse scales. A 
detailed examination of several realizations of such distri- 
butions shows that the optical-depth variations of 0.2-0.4 
across the relatively small tranverse separation (50-100 AU) 
as seen in Fig. 3, are not very rare, but would be expected 
with ~10% probability. 

We would like to re-emphasize here that in opacity dif- 
ference measurements with a transverse separation x , one 
is probing a point of the structure function and not a point 
in the spectrum of opacity distribution. In general, in fact 
all scales in the spectrum contribute to such a measurement, 
except, of course, the uniform component and, ironically, the 
scales x and its integral submultiples that contribute noth- 
ing. In detail, the structure function value at result 
of a sum of the contributions at all spatial frequencies in the 
power spectrum after modulation by 1 — Cos(2irf x x ), and 

considering an ensemble average of such sumsj] correspond- 
ing to all possible orientations and locations of vectors of 
length x ( along with f x , f y axes). This modulation is sim- 
ply a result of the two-point difference measurement. The 
modulating function has its first peak at f x = l/(2x ) and 
at odd multiples of it then on. When the power spectrum is 
red (i.e. a is positive) and steep, even the highly attenuated 
contributions (due to the modulation) from the low spatial 
frequencies (long scales) can, and do indeed, dominate in the 
net contribution, making an equivalent scale being probed 
much larger than x . One may estimate the equivalent scale 
size by considering a weighted average over all scales, where 
weights are determined by the associated values of the mod- 
ulated power spectrum. In a simple estimation, for exam- 
ple, considering a 2-d spectrum over |/ s | < l/x , we find 
the equivalent scale to be nearly one order of magnitude 
longer than x . A proper estimation should include the full 
spectrum. However, since such an equivalent scale has little 
physical meaning, we have not pursued such an estimation 
further. In summary, the earlier interpretations of the opac- 
ity difference observed at a transverse separation X q £LS being 
a result of a cloudlet of the size same as the separation x 
appear, to us, erroneous. 



2.3 The over-dense (& over-pressured) HI 
cloudlets ? 

One of the major mysteries which owes its origin to the 
above mentioned misinterpretation is that of the observed 
opacity changes combined with an assumed longitudinal 
scale implying highly over-dense (uh ~ 10 4 — 10 5 cm -3 ) 
and, consequently, over-pressured cloudlets. Even if the ob- 
served opacity changes were to be accepted at their face 
value, their implication needs to be reinterpreted, since a) 
as emphasized in the previous section, the observed varia- 
tion is contributed by the whole range of scales and not by 
just one particular scale same as the transverse separation 

t More formal analytical expressions involving Bessel functions 
are commonly used. See, for example, Cordes, Weisberg & Bori- 
akoff (1985). 
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probed, and b) the measured value can not be directly asso- 
ciated with a particular longitudinal scale. These considera- 
tions are applicable not just to the "two-point comparisons" 
(i.e. variations expected across a 1-d cut), but also to the 
features observed in the 2-d images of the opacity. The varia- 
tions like those apparent in our 1-d simulation (Fig. 3) would 
be equally probable in an equivalent 2-d image where they 
would appear as one-dimensional features. The edges or the 
elongated features apparent in only some of the 2-d opacity- 
images (Davis, Diamond & Goss 1996; Faison et al. 1998; 
Faison 1999) are therefore not at all surprising, whereas an 
opacity variation feature that is narrow in both dimensions 
should be considered relatively rare. 

If one wants to estimate, using the data, the properties 
of the small-scales in the HI distribution, then the follow- 
ing is one correct way to proceed. One treats the measured 
value of opacity difference or of the associated HI column- 
density change as just an estimate of the associated struc- 
ture function at the probed transverse separation. From this, 
and with some knowledge of the spectrum (or the structure 
function itself), it is possible to estimate the implied power 
in the same scale as the transverse separation probed. Now, 
this power (from the power spectrum) or the amplitude of 
the ripple corresponding to that scale, can justifiably be in- 
terpreted in terms of the associated fluctuating optical depth 
or density. More formally, with the measurement giving an 
estimate of (square-root of) the structure function S Tv (x ), 
the optical depth variation associated with the scale x is 
to be estimated as the (square-root of) power spectral con- 
tribution < P Tv {f s = 7-) >, a value significantly smaller 
than the former. For example, using the relevant values ob- 
served in the Cas-A direction, a observed change of 0.2 in 
the optical depth between two sight-lines with a transverse 
separation of about 1000 AU, would imply a r fluctuation 
on the scale of ~1000 AU to be about 10~ 5 . Assuming a 
velocity-channel width of 0.5 km/s & a spin temperature of 
about 100 K, the contribution (or deficit) in the volume- 
density from that scale would be less than 0.1 cm -3 , very 
much smaller than what the earlier interpretations would 
suggest ! This value would be even smaller when the possi- 
ble statistical enhancement due to the finite thickness of the 
medium is accounted for. For example, if the contribution is 
from more than one, say, N layers along the sight-line, then 
the corresponding contribution to volume density would be 
y/~N times smaller. It follows that the well understood parsec 
scale would contribute an HI volume density of ~ 1 cmT 3 
or smaller, in good agreement with relevant observations. 
Of course, the actual density at given spatial location would 
be a sum total of such contributions also from a hierarchy 
of scales longer and shorter than x . It is easy to show that 
such contributions to volume density at a given spatial point 
would follow a power law as a function of the scale-size. For 
example, in the Cas-A direction, it may be expressible as 

Ang(i ) ~ Ax a 2 , where x is expressed in parsec, A is a 
constant close to unity and (3 is between 1 to 2 depending 
on whether our sight-line encounters contributions many (of 
the order of l/x , with x in pc) or just one layer of scale x , 
respectively. Even in the worst case, i.e. when /3 = 2, the rms 
fluctuations in the volume density, estimated on the parsec 
scales, would be below 100 cm -3 . The more likely value of (3 
in the sub-parsec regime of scales is closer to unity, and then 



the volume density fluctuations would be moderate, with an 
rms of ~10 cm -3 . The detailed quantitative picture may 
differ between different directions in the Galaxy by a fac- 
tor of 10 or less, and would hopefully be revealed by future 
suitable observations & a careful interpretation. 

To conclude, there appears to be no compelling obser- 
vational evidence for the so called "highly over-dense small 
scale structure", and even the observations probing small 
transverse scales are not at all inconsistent with what we 
would expect by extrapolating from the better studied range 
of large and moderate spatial scales in the cold neutral 
medium. 

2.4 The reported measurements versus 
uncertainties 

So far, we have taken the reported observations (estima- 
tions) of the opacity differences at their face value. How- 
ever, certain uncertainties inherent to the measurements are 
worth noting. As we have discussed earlier, any observed 
line-visibility in VLBI observation results from structure in 
both the background source and the opacity as seen by a 
given interferometer baseline. The interpretation can be- 
come more complicated when only a limited number of base- 
lines are used and can even be misleading if any changes 
in the relative orientation of the baselines (in addition to 
its projected length) are not accounted for. As was already 
pointed out (Radhakrishnan & Deshpande 1990; unpub- 
lished), the uniqueness of the interpretation of Diamond et 
al. (1989) becomes debatable on these grounds. Later sim- 
ilar observations however have resolved the possible ambi- 
guities by actually mapping the opacity distribution across 
the extra-galactic sources and hence the estimated opacity 
differences can be considered reliable. 

One general but important aspect, addressed earlier by 
Deshpande et al. (1992), is that of the contribution from 
the HI emission to the measurement uncertainty (relevant 
to both the interferometric and pulsar probes). Even if a 
VLBI baseline resolves out the large scale HI emission, each 
of the elements of the interferometer does respond to the HI 
emission contribution and the equivalent system tempera- 
ture of the interferometer can be significantly higher in the 
corresponding spectral channels compared to that for only 
continuum contribution. As for the single dish observation, 
the HI emission contribution to the estimation uncertainty 
is rather obvious. The measurements in the pulsar directions 
(Deshpande et al. 1992; Frail et al. 1994) as well as more re- 
cent interferometric measurements (Davis, Diamond & Goss 
1996; Faison et al. 1998; Faison 1999) do explicitly take into 
account the HI emission contribution. Interestingly, the sys- 
tem temperatures in the line channels also depend, in prin- 
ciple, on the optical depths at the corresponding frequen- 
cies. This can be a significant effect when the background 
sources make dominant contribution to the system temper- 
ature (in continuum channels). In such cases, higher optical 
depths result in significantly smaller system temperature in 
the corresponding spectral channels. 
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3 DISCUSSION 

In this paper we have addressed some aspects related to 
the observation, analysis and interpretation of opacity differ- 
ences across small (sub-parsec) transverse scales. Although 
the considerations we raise are simple-minded, they appear 
to have serious implications that argue against certain inter- 
pretations, such as those suggesting the so called small-scale 
structure of highly over-dense cold HI cloudlets as being re- 
sponsible for the observed opacity differences. 

We have emphasized the need for recognizing the na- 
ture of the actual quantity one measures through the probes 
that have been employed and that almost all scales con- 
tribute to the measured opacity differences. But, so far, the 
existing studies have misinterpreted the observed opacity 
difference between two sight-lines of the associated trans- 
verse separation as due to opacity structure on that scale 
in three-dimensions, and therefore, it is not surprising that 
the implied volume densities appear extra-ordinarily large. 
As illustrated in the earlier sections of this paper, the ob- 
servations appear consistent (both, qualitatively and quan- 
titatively) with a single power-law description of the HI dis- 
tribution over the entire relevant range of scales and do not 
imply any mysterious structure. In our simple picture, the 
spectral behaviour studied up to moderate scales (e.g. ~ 0.1 
pc scale as by DDG, 2000) is assumed to extend with the 
same power-law index to the 10 AU scale. We are aware of 
the study by Croviser, Dickey & Kazes (1985) that claimed 
a cutoff in structures below 0.2 pc. Considering 1) their 
method in which signatures of small-scale structure were 
expected to show up close to the zero- velocity, 2) the major 
absorption line features (differences in which were probed) 
were well away from zero-velocity, and 3) the velocity res- 
olution was coarse, we think that their study suffers from 
serious selection against structures smaller than about 0.5 
pc. 

The recent suggestion by Heiles (1997) did for the first 
time distinguish between the transverse scale probed and the 
longitudinal scale for estimating the volume density, but by 
invoking "thin" structures such as filaments & sheets that 
should preferentially align along sight-lines. As we have dis- 
cussed, the distinction between the transverse scale probed 
and the corresponding equivalent longitudinal scale appears 
to be rather simple and more inherent to the basic measure- 
ment than any geometrical shapes would imply. While any 
anisotropy in structures is likely to increase (statistically) 
the expected magnitude of the opacity differences (and elon- 
gation factors up to ~2 may be common) it does not appear 
necessary for understanding the available observations. In 
any case, the structures suggested by Heiles may be difficult 
to produce and maintain, particularly in the certain align- 
ment that they need to have with our sight-lines. 

The velocity spread associated with turbulence can 
produce additional corrugations in the opacity distribution 
when viewed over velocity channels narrow compared to the 
spread due to turbulence (as is normally the case). It would 
be instructive to assess such an effect quantitatively. 

It may be important to note that the radio/optical ob- 
servations of opacity changes in other species (e.g., a. H2CO 
and OH by Moore & Marscher 1995; b. Nal and Call by 
Mayer & Blades 1996, and Watson & Meyer 1996) should be 
interpreted much the same way as we have discussed in the 



context of HI and then a similar puzzle these observations 
appeared to have raised should stand resolved. The discus- 
sion in this paper is relevant also to the dispersion measure 
changes detected in pulsar directions and in general to any 
situation involving a similar probe. 

Although in the course of our discussion we have drawn 
upon the HI data in the Cas-A direction as a useful example, 
we do recognize the possibility that the atomic medium in 
other directions may have quite different column densities as 
well as power-spectral descriptions from that in the Cas-A 
direction. However, the main issues we have addressed are 
of a more general nature and are not crucially based on the 
quantitative estimates from the Cas-A data. The HI emis- 
sion line study by Green (1993) indicates that the power 
spectra (and structure functions) may be less steep in some 
directions than that for cold HI in Cas-A direction. A recent 
study of HI in the Small Magellanic Cloud suggests a rela- 
tively steeper power spectrum (Stanimirovic et al. 1999). If 
the power spectra of scale distribution in cold HI also have 
a similar variation, then in some directions we should ex- 
pect even larger opacity differences on small transverse sep- 
arations (much more than even those seen in recent VLBI 
studies). For example, a change of 0.1 in a would increase 
the expected opacity changes on the AU scales by a factor of 
~2. 

Further investigations should benefit from using the 
available data for a systematic estimation of the structure 
function associated with the opacity distribution over the 
relevant range of transverse scales. HI absorption measure- 
ments on moderate size, bright background sources should 
help extending the direct power spectral analysis (e.g in the 
case of Cas-A) to intermediate and small transverse scales. 
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